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PROCEDURE 
Calibration Tests 

A number of preliminary investigations were med* in order to 
check the Instrumentation. 

Air -flow measurements . - The flow nozzle at the heater-tube 
entrance was calibrated for the unheated condition by using the 
large surge tank as an air reservoir. Air was allowed to flow from 
the surge tank through the nozzle; the pressure at the nozzle 
entrance was held constant. The pressure temperature in the 
surge tank were recorded at selected time intervals; by use of the 
equation of state, the flow rate was obtained. The flow rates so 
obtained were from 0 to 2 percent higher than the values calculated 
from A.S.M.E. equations. The rotameter, which had been calibrated 
by a gasometer, showed about the same deviation from the calculat ed 
flows. 


Tube -wall -temperature measurements . - In order to determine 
whether errors existed in tube -wall -temperature measurements due 
to the current in the heater tube, power was supplied to the tube 
with no air flowing until a predetermined wall temperature, as 
indicated by the thermocouples, was obtained. The power was then 
shut off and the readings of several thermocouples were recorded at 
5-second intervals . The temperature s determined by extrapolation to 
zero time of the temperature against time curve thus obtained for 
each thermocouple agreed with the temperatures read Just prior to the 
power shutoff; hence it was concluded that the current had no effect 
on wall -temperature measurements. 

Heat losses. - The heat loss from the tube was obtained for the 
range of tube -wall temperature by supplying various amounts of 
power to the heater tube with no air flowing through the system. 

After equilibrium conditions had been maintained for approximately 
l/2 hour, the power input and all tube -wall temperatures were 
recorded. The power input for a given average wall temperature with 
no flow was considered to be the heat loss for the same average wall 
temperature with air flowing in the tube. 


Heat-Transfer Tests 

investigations were conducted to obtain surface heat -transfer 
coefficients and the associated static-pressure drops for average 
inside -tube -wall temperatures from 220° to 1240° F (corresponding 
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maximum local wall temperatures, 250° to 1500° F), Reynolds numb ers 
fron 10,000 to 250,000, and the tube -exit Mach numbers up to 1.0. 

The air temperature at the tube entrance was about 85° F and the 
pressure varied from 15 to 45 pounds per square inch absolute. 

The experimental procedure was as follows: The entrance -air 
pressure In the calming tank was set at the minimum value and power 
was supplied to the heater tube until the desired tube -wall tempera- 
ture was obtained. After equilibrium conditions had been maintained 
for approximately l/2 hour, all power input, pressure, and temper- 
ature readings were recorded with the Franz thermocouple probe 
retracted in the outlet mixing tank. The probe was then moved 
upstream along the center line of the tube and the air temperature 
was recorded at 4-inch intervals . This procedure was repeated 
over the range of entrance pressure for five wall-temperature levels. 


SYMBOLS 

The following symbols are used in the calculations: 
c p specific heat of air at constant pressure, (Btu/(lb)(°F)) 

1) inside diameter of heater tube, (ft) 

G mass velocity of air, (lb/(hr)(sq ft)) 

h heat-transfer coefficient, (Btu/(hr)(sq ft)(°F)) 

k thermal conductivity of air, (Btu/(hr)(sq ft)(°F/ft)) 

kj thermal conductivity of Inconel, (Btu/(hr)(sq ft)(°F/ft)) 

L effective heat -transfer length of heater tube, (ft) 

dp static -pressure drop across heater tube, (lb/sq ft) 

Q rate of heat transfer to air, (Btu/hr) 

Ri inner radius of heater tube, (ft) 

R 0 outer radius of heater tube, (ft) 

S heat -transfer area of heater tube, (0.211 sq ft) 

Tq total air temperature at entrance to heater tube, (°F) 
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total air temperature at exit of heater tube, (°F) 
average bulk temperature of air, ^(T^ + T 2 ), (°F) 
average Inside -vail temperature of heater tube, (°F) 
average out side -vail temperature of heater tube, (°F) 
bulk velocity of air, (ft/hr) 
air flow, (lb/hr) 

absolute viscosity of air, (lb/(ft)(hr)) 
density of air, (lb/cu ft) 

ratio of densities at entrance and exit of heater tube 

ratio of density at tube entrance to standard sea-level 
density, p-j/0.0765 

Prandtl number 

Nusselt number 

Reynolds number 

P 0 V t> D/p s modified Reynolds number 

Subscripts : 

b physical properties of air evaluated at average bulk 

temperature (average total temperature T^) 

s physical properties of air evaluated at average inside- 

wall temperature T s 

METHODS OF CALCULATION 

Temperatures , - The average out side -tube -vail temperature T 0 
was obtained by plotting curves of -temperature against tube length, 
measuring the area under the curve, and dividing by the tube length. 
The plotted values of temperature were the average of the two thermo 
couples located 180° apart at each of the 15 stations along the tube 


‘s 


W 

P 

P 

Pl/P2 

a l 


Cpp/k 

hD/k 

DG/p 
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The average Inside -tube -vail temperature T e was calculated, 
from the outside -tube -wall temperature, the heat flow, and. thermal 
c onduc tivity and physical dimensions of the tube by the following 
equation, which can be derived (reference 1) with the assumptions 
that heat is uniformly generated across the tube wall and that heat 
flow is directed radially inward: 


T s “ T o 


- Bi 2 ) 



V B o 

log eRr- 



When the dimensions of the heater tube and a value of 6.6 
for the thermal conductivity of Inconel are substituted in this 
equation 


T b = T o - 0.00107 Q 

A constant value of thermal conductivity was used because of 
lack of data on variation with temperature. The error introduced 
is negligible Inasmuch as the difference between inside- and 
outside-wall temperature is small compared vith the difference 
between average wall and air temperature (for example, maximum 
value - 23° F temperature difference between outer and inner wall 
at 940° F temperature difference betveen wall and air) . 

The average bulk air temperature Tj, was taken as the arith- 
metic mean of the total temperatures and T 2 of the air meas 
ured in the inlet and the outlet tanks. 

Heat-transfer coefficients. - Average heat-transfer coeffi- 
cients h were calculated from the following equation: 


h = 



S 



(t 2 - V 

•h * Hi ) 


W Vbt T 2 - V 

s 1 t s - v 


The use of average air temperature in computing the heat- 
transfer coefficient is considered satisfactory inasmuch as the 
temperature rise of the tube wall along its length; neglecting 
the end points, was about equal to the rise in air temperature 
from inlet to outlet, as will be shown later. 
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Physical properties of air. - The specific heat at constant 
pressure c p , the thermal conductivity k, the viscosity u, and 
the corresponding Prandtl number Cp|i/k used were taken from ref- 
erence 2. These physical properties are plotted in figure 4 as 
functions of temperature. 


RESULTS AND DISCUSSION 
Temperature Distribution 

The axial distribution of outside -tube -wall temperatures along 
the tube length for the condition of no air flow is shown in fig- 
ure 5 for various amounts of electrical heat input. High tempera- 
ture gradients exist at the entrance and the exit of the tube as a 
result of end conduction losses through the electrical connector 
flanges and cables; the temperature, however, is reasonably constant 
along the greatest portion of the tube length. The slight dips, 
which occur in the curves between 8 and 12 inches from the tube 
entrance, are not understood at this time. 

Representative axial distributions of out side -tube -wall tem- 
peratures are shown in figure 6 for five different amounts of total 
electrical heat input at an approximately constant average Reynolds 
number of 140,000. The Reynolds number was computed at the average 
air temperature T^, hence the air flow had to be increased with 
increasing heat input and attendant increasing exit and average air 
temperature in order to maintain constant Reynolds number. The 
temperature increases almost linearly along the central 20 inches 
of the tube and, similar to the no-flow case (fig. 5), decreases 
sharply at both ends of the tube. 

The distribution of unraixed air temperature along the center 
line of the tube, as indicated by the Franz thermocouple probe, is 
shown in figure 7 for the same conditions as figure 6. Included 
in the figure at the 24-inch station are the mixed exit-air temper- 
atures as measured in the mixing tank. The unmixed air temperature 
increases almost linearly along the central 16 to 20 inches of tube 
length and levels off at both ends because of decreased tube -wall 
temperature (fig. 6). The mixed exit-air temperature is higher 
than the unraixed air temperature indicated by the probe at the 
24-inch station; the difference in temperature increases from 17° F 
at the lowest heat input to 180° F at the highest heat input. These 
differences in temperature might be explained by the existence of 
large temperature gradients across the air stream in the tube. 
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Examination of figures 6 and 7 indicates that the rise in mixed 
air temperature (fig. 7) is approximately equal to the increase in 
wall temperature along the tube length if end effects are neglected 
(fig. 6). 


Heat Balance 

The external heat losses from the heater tube as obtained from 
figure 5 are cross-plotted against average outside-wall temperature 
of the heater tube in figure 8. The average outside-tube-wall tem- 
peratures were obtained by integration of the curves of figure 5. 

A heat balance for all the heat -transfer experiments is shown in 
figure 9 in which values of the heat transferred to the air 

[Wc p ,b(T2 - Tx)j plus the heat lose (obtained from fig. 8) are 
plotted against electrical heat input. 

A good heat balance is indicated with a maximum deviation of 
about 5 percent from the match line. 


Correlation of Heat-Transfer Coefficients 


Correlation based on bulk temperatures. - Forced convection, 
turbulent flow, heat-transfer coefficients for low viscosity 
fluids are generally correlated by means of the familiar Nusselt 
relation (reference 3) where Nusselt number divided by Prandtl 


number to a power (generally 0.4) M 





0.4 


is plotted 


against Reynolds number (DG/n^) 8114 in which the physical prop- 

erties of the fluid are evaluated at the average fluid bulk tem- 
perature. The results of the present investigation are plotted in 
this manner in figure 10. Included for comparison is the average 
line (dashed) obtained in reference 3 from a correlation of the 
results of various investigators and the average line (dotted) 
corresponding to the correlation obtained in reference 4 for the 
data of the same investigators using the same physical properties 
(reference 2) used herein (fig. 4). The dotted line was obtained 
by transposition of the Stanton number plot of reference 4 to the 
coordinates of figure 10 wherein Prandtl number was assigned a 
value of 0.70 corresponding to the physical properties of air (fig. 
at a temperature of 120° F. The line corresponding to the corre- 
lation of reference 4, although obtained for the same data, is 


4 ) 
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lower than that of reference 3, undoubtedly because of use of differ- 
ent physical properties. The lower line is considered preferable for 
comparison with the correlation of the data of this investigation 
inasmuch as the same physical properties are used. 


A family of parallel lines for the different temperature levels 
is obtained for the present data with slopes of approximately 0,8 
at Reynolds numbers above approximately 25,000. At lower Reynolds 
numbers the data fall off indicating the'presence of transition. 

The low temperature data (220° F wall temperature) are in good agree- 
ment with the line of reference 4; however, the higher temperature - 
data lines decrease progressively with increased wall temperature and 
corresponding increased temperature difference between wall and air. 


The values of 



for constant Reynolds number 


decrease about 25 percent with increase in average wall temperature 
from 220° to 1240° F. 


Correlation based on surface temperature. - The data of fig- 
ure 10 are replotted in figure 11 wherein the physical properties 
are evaluated at the inside -tube -wall or surface temperature instead 
of the air bulk temperature. This plot shows a similar separation 
of the data with temperature level as obtained in figure 10 except 


that the separation is increased; the values of 



decrease approximately 40 percent as the tube temperature level 
increases from 220° to 1240° F. If the physical properties were 
evaluated at the film temperature (average of the surface and bulk 
temperatures), the separation of the data with temperature level 
would be between the 25 and 40 percent indicated in figures 10 and 
11, respectively. 


Correlation based on modified Reynolds number. - A plot of 
the data is shown in figure 12 wherein the physical properties are 
again evaluated at the surface temperature but in addition Reynolds 
number is modified by substituting' the product of air density 
evaluated at the surface temperature p s and average air velocity 
evaluated at the air bulk temperature for the conventional 

mass velocity G or In addition, the viscosity term in 

Reynolds number is evaluated at the surface temperature . The rela- 
tion between the modified and the conventional Reynolds number is 
given by the following equation: 
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p s Y b D _ /dg\/£s\ /Do V^b * 460> \ _ foq\ Pb * 46Q \ fa>\ 

H s " WvV = VsA T b + 460 / = W\ T s + 460 / Vs/ 

This definition of Reynolds number is essentially similar to that 
theoretically shown in reference 5 to be the correct value for 
laminar heat transfer from a flat plate. This method of plotting 
(fig. 12) results in a good correlation of the data for Reynolds 
numbers above transition for all the temperature levels Investigated 
(average wall temperatures up to 1240° F and corresponding maximum 
local wall temperatures up to 1500° F) . The equation of the line, 
which represents the data for Reynolds numbers above approximately 
10,000, is 



0.022 


Pfl^bD 
\ H e ' 


, 0.8 


In figure 13 the data are re plotted neglecting Prandtl number, 
that is, hD/k s is plotted against p- . A good correlation 

is again obtained indicating that the variation of Prandtl number 
vith temperature (fig. 4) has little effect. The equation of the 
line best representing the data is 


££ = 0.018 
o 


/p s V D 


, 0.8 


Pressure -Drop Correlation 

The measured static-pressure drops across the heater tube are 
correlated in figure 14 where the product of pressure drop and 
entrance density (referred to standard sea-level density) divided 

0 1 Ap 

by air flow to the 1.8 power y •g 1 iB plotted against the ratio 

w 1.8 

of entrance -to-exit density P 1 /P 2 ♦ The exponent of 1.8 on air 

flow was obtained from preliminary logarithmic plots of OjAp 
against W for constant values of p^/pg. The air densities p^ 
(and corresponding 0^) and Pg were calculated using static 
pressure and total temperature at entrance and exit of the tube, 
respectively. Total instead of static temperature was used in 
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order to simplify the calculations involved in reducing the experi- 
mental data. All the pressure-drop data correlate within a devia- 
tion of about ±1 percent except for four points, which were obtained 
at Beynolds numbers below 10,000. The equation of the line through 
the data is 

= 0.10 - 0.044 

W 1 * 8 p 2 


SUMMAEY OF RESULTS 

The results of the heat-transfer investigation conducted with 
air flowing through an electrically heated Inconel tube with a 
rounded entrance, an Inside diameter of 0.402 inch, and a length 
of 24 inches over a range of Reynolds numbers up to 250,000 and 
a range of average inside -tube -wall temperatures up to 1240° F 
(corresponding maximum local tube -wall temperatures up to 1500° F) 
showed that: 

1. Correlation of the average heat-transfer coefficients 
according to the familiar Nusselt relation wherein physical proper- 
ties were evaluated at the average bulk temperature resulted in a 
•separation with temperature level of the parameter Nusselt number 
divided by Prandtl number to the 0.4 power; at constant Reynolds 
number, the parameter decreased about 25 percent with increase in 
average surface temperature from 220° to 1240° F. Evaluation of 
the physical properties at the average surface temperature reduced 
the parameter by 40 percent for the same increase in temperature. 

2. A good correlation of the heat -transfer data was obtained 
for the entire range of temperatures when the viscosity |i B , 
thermal conductivity k s , and the specific heat Cp^ s of the air 

were evaluated at the average surface temperature and the Reynolds 
number was modified by substituting the product of air density 
evaluated at the average surface temperature p B and velocity 

evaluated at the average air bulk temperature Vjj for the conven- 
tional mass velocity G. The equation of the line representing 
the data for Reynolds numbers above 10,000 is 
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Satisfactory correlation vas also obtained neglecting Prandtl num- 
ber in which case the equation for the line fitting the data is 


hD 

k s 


0.018 




v0.8 


^8 


3. The pressure-drop data were correlated by plotting the 
product of static -pressure drop Ap a n d entrance air density 
(referred to standard sea-level density) ,0]_ divided by air flow 
W raised to an exponent against the ratio of entrance -to-exit 
density p-^/pg • The equation representing the data within a 

±1 percent deviation for Reynolds numbers over 10,000 is 

^lAp P x 

Ta = 0.10 T? - 0.044 

w 1.8 P2 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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Figure 2. - Schematic diagram of test section showing thermocouple and pressure-tap locations 
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Figure 3. - General view of heater-tube installation 
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Figure 4. - Variation of specific heat Cp, thermal conductivity k, absolute viscosity U, 
and Prandtl number Cpu/k of air with temperature. (Data from reference 2.) 
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Figure 5. - Axial distribution of outside-tube-wall temperature for 
various amounts of electrical heat input with no air flow* 
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Figure 6. - Representative outside-tuba-wall temperature distributions 
for various amounts of electrical heat input and a constant Reynolds 
number of approximately 140,000. 











External heat loss, 6tu/hr 


NACA RM NO. E7L3I 


23 



Figure 8. - Variation of external heat loss with average outside-tube- 

wall temperature. 


Heat transferred to air + external heat loss, Btu/ 


24 


NACA RM NO.' E7 L 3 1 



Figure 9 


Heat balance 










500 















HACA EM No. E7L31 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

HEAT TRANSFER FROM HIGH -TEMPERATURE SURFACES TO FLUIDS 
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WITH ROUNDED ENTRANCE, INSIDE DIAMETER OF *0.4 INCH 

AND LENGTH OF 24 INCHES 
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*\s , - 

SUMMARY 

-..y 

A preliminary heat -transfer investigation was conducted vith 
air flowing through an electrically heated Inconel tube vith a 
rounded entrance, an inside diameter of 0.402 inch, and a length of 
24 inches over, a range of Reynolds numbers up to 250,000 and a 
range of average inside -tube -wall temperatures up to 1240° F. The 
corresponding meu Imum local tube -wall temperatures ranged up to 
1500° F. 

Correlation of the heat -transfer data by the conventional 
method wherein physical properties of the air were evaluated at the 
average bulk temperature resulted in a reduction of Nusselt number 
of about 25 percent for an increase in average surface temperature 
from 220° to 1240° F at constant Reynolds number. A good correla- 
tion of the data for the entire temperature range was obtained, 
however, when the physical properties of the air were evaluated at 
the average surface temperature and the Reynolds number was modi- 
fied by substituting the product of air density evaluated at the 
average surface temperature and velocity evaluated at the average 
air bulk temperature for the conventional mass flow per unit cross- 
sectional area. 

Static pressure drops were obtained for all operating condi- 
tions and were satisfactorily correlated with air flow rate and 
the ratio of entrance -to -exit density. 
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INTRODUCTION 

Moat of the fundamental data available on forced convection 
heat transfer between surfaces and fluids have been obtained at 
relatively low temperatures a n d heat-flux densities and do not 
extend into the range of high temperatures and fluxes that are of 
interest in many current engineering applications. The accuracy of 
extrapolation of the existing data to high-temperature and high -flux 
conditions is doubtful because of the possible existence of extreme 
temperature and velocity gradients in the fluid film, changes in 
physical properties of the fluid, and radiation effects that may 
became important for the case of surface -to-gas heat transfer; 
changes in boiling phenomenon will markedly affect surface -to- 
liquid heat transfer. 

Accordingly, an experimental investigation has been instituted 
at the NACA Cleveland laboratory to obtain surface -to-fluid heat 
transfer and associated pressure -drop information over a range of 
temperatures and heat fluxes. As part of the general program, an 
investigation is being conducted with air flawing through an elecr 
trically heated Inconel tube having an inside diameter of 0.402 inch 
and an effective heat -transfer length of 24 inches.- Experiments 
have been conducted with a rounded tube entrance over ranges of 
Reynolds numbers up to 250,000, average tube -wall temperatures up to 
1240° F with corresponding maximum local tube-wall temperatures up 
to 1500° F, heat fluxes up- to 100,000 Btu per hour per square foot, 
and tube -exit Mach numbers up to 1.0. 

The results of these preliminary investigations, which were 
obtained during the fall of 1947, are reported herein. Average heat- 
transfer coefficients are correlated in accordance with the familiar 
Nusselt relation and with a modification thereof; pressure drops are 
correlated with air-flow rate and density changes. The correlation 
of these preliminary data, although very satisfactory, will be sub- 
ject to further investigation over a greater range of conditions. 


APPARATUS 

A schematic diagram of the heater tube and the associated 
equipment used in this Investigation is shown in figure 1. 


Heater Tube 

The details of the heater tube are shown in figure 2. The 
test section consists of an Inconel tube having an inside diameter 


